High energy photons from passage of jets through quark gluon plasma 
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We calculate the production of high energy photons from Compton scattering and annihilation of 
a quark jet passing through a quark gluon plasma produced in a relativistic heavy ion collision. The 
contributions are large and reflect the momentum distribution of the jets and the initial conditions 
of the plasma. 



Relativistic heavy ion collisions are studied with the 
aim of producing a plasma of quarks and gluons (QGP). 
Photons are considered to be an important probe for 
the investigation of the formation and evolution of such 
a plasma due to their weak final-state interactions Q. 
Once produced they carry the information about the con- 
ditions of the environment in which they were created, 
encoded in their momentum distribution, thus providing 
a glimpse deep into the bulk of strongly interacting mat- 
ter. Though most of the measured photons have their 
origin in the decay of hadrons after the QGP phase, it 
has become possible to isolate the direct photons pro- 
duced in such collisions 0. 

The sources of direct photons considered so far 
include quark annihilation, Compton scattering, and 
bremsstrahlung following the initial hard scattering of 
partons of the nuclei jSl, as well as thermal photons from 
the QGP 0, 0, 0, Slal and from hadronic interactions 
in the hot hadronic gas after the hadronization of the 
plasma 0, 0|. The pre-equilibrium production of pho- 
tons has also been investigated by several authors [lOj. 
Results are available for production from the entire his- 
tory of the system ^ . 

In this letter we study a new source of direct photons 
originating from the passage of the produced high energy 
quark jets through the QGP (jet-photon conversion). A 
fast quark passing through the plasma will produce pho- 
tons by Compton scattering with the thermal gluons and 
annihilation with the thermal antiquarks. This process 
is higher order in as compared with photons from initial 
hard scatterings, but it is not a subleading contribution, 
since it corresponds to double scattering, which is en- 
hanced by the size of the system. For cold nuclear mat- 
ter this effect is encoded in multi-parton matrix elements 
which are enhanced by powers of A^^^ ^^l- Below, we 
find that this source is at least comparable in strength 
to the other direct photon sources and even dominates 
in the range p± < 6 GeV for Au-|-Au collisions at the 
Relativistic Heavy Ion Collider (RHIC). 

We also demonstrate that the p_L-distribution of these 
photons is directly proportional to the momentum dis- 



tribution of jets at an early stage after their production, 
before they have lost energy on their travel through the 
plasma. Since the measured high-pj^ hadron spectrum is 
proportional to the spectrum of partons after they have 
left the plasma, a comparision of both spectra could pro- 
vide a quantitative determination of the ener gy l oss and 
help confirm the mechanism of jet quenching \l?L . 

Furthermore, the photon yield depends on the inte- 
grated density of the matter traversed by the jets and 
thus can provide a measurement of this quantity. We 
emphasize that our mechanism is distinct from the back- 
to-back correlation among direct photons and leading 
hadrons, which was proposed as a direct measurement 
of the jet energy loss in dense matter 15]. 

The idea of jet-photon conversion is based on the prop- 
erties of the annihilation and Compton cross sections. 
The kinematics of the annihilation of a quark antiquark 
pair {q-\-q—^ "f + g) is expressed in terms of the Mandel- 
stam variables s — {pq-\-pq)'^, t — {pq—p^Y and u — {pq — 
Pj)'^, where the pi are the four-momenta of the particles. 
The differential cross section at Born level for massless 
partons has the form da/dt = SiraasC^ (u/t + t/u) /9s^, 
where Cq is the fractional charge of the annihilating 
quarks. This implies that the largest contribution to the 
production of photons arises from small values of t or 



u, corresponding to p 



Pq or p. 
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Pq. The pro- 



cess can be visualized as a conversion of one of the an- 
nihilating quarks to a photon. This picture immediately 
suggests that the photons from this process provide a 
direct measurement of the quark (or antiquark) momen- 
tum. For the Compton process {g + q ^ + q) we have 
da/dt — — Tracts + s/'^J-) /3s^ . The dominant con- 
tribution now comes from the region of small u, when 
P7 ~ Pqi corresponding to the conversion of the quark 
(or the antiquark) into a photon. 

Following 16,LL3| we approximate the invariant photon 
differential cross sections for the annihilation and Comp- 
ton processes as 



da^< 
d^p^ 



s)]^E^[5{p^ 



Pq) +'5(P7 -Pq)] (1) 
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i?7^5^«'^^ Hs)i?7'5(P7-Pq)- (2) 

Herecr('^Hs) and cr(c)(s) are the corresponding total cross 
sections. Now consider an ensemble of fast quarks passing 
through a hot medium. A certain fraction, determined 



by the total cross section, will undergo these annihila- 
tion or Compton processes initiated by the medium and 
will emit photons. Whenever this happens the energy of 
the emitted photon will reflect the initial energy of the 
converted quark. 

Using I|1I2|I the rate of production of photons due to 
annihilation and Compton scattering are given by 



'¥x¥p^ ~ 2(2 



Mf 



q- 

Nf 



d^x (Pp-, (27r)6 ^ 
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The fi are distribution functions for the quarks, an- 
tiquarks, and gluons. Inserting thermal distributions 
for the gluons and quarks one can obtain an analytical 
expression for these emission rates for an equilibrated 
medium 0,0, E3- 

We propose that the phase-space distribution of the 
quarks and gluons produced in a nuclear collision can be 
approximately decomposed into two components, a ther- 
mal component fth characterized by a temperature T and 
a hard component /jot given by hard scattering of par- 
tons and limited to transverse momenta p± ^ 1 GeV/c: 
/(p) = /th(p) + /jot(p)- We note that /jct dominates for 
large momenta, while at small momenta / is completely 
given by the thermal part. 

The phase-space distribution for the quark jets prop- 
agating through the QGP is given by the perturbative 
QCD result for the jet yield (191] : 



/jot(p) 



1 (2^) 



JCt 



gq ttR\tp± cPp±dy 

X (5(77 - - Ti)e(T,. 



Rir) 
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where gq — 2 x 3 is the spin and colour degeneracy of 
the quarks, R± is the transverse dimension of the sys- 
tem, Ti l/p± is the formation time for the jet and t] is 
the space time rapidity. R{r) is a transverse profile func- 
tion. We shall take Tmax as the smaller of the life-time Tf 
of the QGP and the time taken by the jet produced 
at position r to reach the surface of the plasma. The 
boost invariant correlation between the rapidity y and 77 
is assumed here 0| . Similar results are obtained when a 
Gaussian correlation between y and ri, characterized by 
a width Ap ~ 2/p± coshy is assumed [l9l| . 

We calculate the jet production and the direct pro- 
duction of photons in lowest order pQCD We have 
used the CTEQ5L parton distributions [20] and EKS98 



nuclear modifications [2l| with the scale set to p±. We 
further use a X-factor of 2.5 in the jet production to ac- 
count for higher order corrections, cf. 22]. We neglect 
the weak dependence of K on p± and \/S in this study 
1 2 3*] . Our results for the number of quarks, antiquarks 
and gluons at y = can be parameterized as 



dm 



jot 



d'^p± dy 



= T, 



du^" 



AA 



y=0 



d'^pj_ dy 



Here Taa — 9A'^ /SttR'^ is the nuclear thickness for a 
head-on collision. Numerical values for the parameters 
a, b and c are listed in Tabled 

The integral over p in Eqs. Q and is dominated 
by small momenta. Wc can therefore drop the jet part in 
the decomposition of the distributions /(p) in the inte- 
grands and approximate them by the thermal part. After 
performing the integrals the results for annihilation and 
Compton scattering are identical and read [l^ 
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with C = —1.916. We have included the three lightest 
quark flavours so that J2f ^g//^^ ~ 2/3. 

Once again we apply the decomposition of the distri- 
butions to 10. Inserting the thermal part we recover 
the known rate of emission of high energy photons from 
interactions within the plasma. Since we are dealing 
with (nearly) massless partons we have to address the 
infrared divergence in (jZj) in the limit m — > 0. It has 
been shown jifS] that the divergence can be eliminated 
by including higher-order effects from the interaction of 
the quarks in the thermal medium with the result that 
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a [l/GeV^] 


b [GeV] 


c 


RHIC 


q 


5.0x10^ 


1.6 


7.9 


q 


1.3x10^ 


1.9 


8.9 


LHC 


q 


1.4x10'' 


0.61 


5.3 


q 


1.4x10'^ 


0.32 


5.2 



TABLE I: Parameters for the minijet distribution dN/cPp± dy 
eiven in Eq. at 1/ = f or Au +Au at ^ft^ = 200 GeV 
(RHIC) and for Pb+Pb at y^S^ = 5.5 TeV (LHC). Numbers 
for quarks and antiquarks are mean values for the three light- 
est flavours: q=(u-|-d+s)/3, q=(u+d+ s)/3. K = 2.5. The 
range of vahdity is 2 GeV/c < px < 20 GeV/c. 



TO in the argument of the logarithm must be replaced 
with 2tojj^ = g^T^ /3. We assume that these considera- 
tions remain valid for the emission of photons even when 
one of the partons is from a jet while the other one is 
from a thermal medium. This assumption remains to 
be verified. We add that the thermal photon production 
from the plasma has recently been calculated up to two 
loops fi\ and to complete leading order 0, and is about 
a factor of two larger than the lowest order results given 
here. It will be of considerable interest to extend our 
present work to account for these effects. 

Now the rate of high energy photon emission from the 
passage of quark jets through the QGP is obtained by 
substituting the jet contributions for fq and /g in 0. In 
order to calculate the photon yields we specify our model 
for the medium. Here we are interested in the emission 
of photons from the QGP phase alone. We assume that 
a thermally and chemically equilibrated plasma is pro- 
duced in the collision at time tq with temperature To, 
and we ignore the transverse expansion of the plasma. 
Further assuming an isentropic longitudinal expansion 
10, we can relate TqTq to the observed particle rapid- 
ity density dN/dy. We take dN/dy « 1260, based on 
the charged particle pseudorapidity density measured by 
the PHOBOS experiment 2A] for central collisions of 
Au nuclei at ^/Snn = 200 GeV. For central collision 
of Pb nuclei at LHC energies we use dN/dy = 5625 es- 
timated by [2^. Imposing a rapid thermalization lim- 
ited by To ^ l/3To |25j we fix our initial conditions to 
be To = 446 MeV and tq = 0.147 fm/c for RHIC and 
To = 897 MeV and To = 0.073 fm/c for LHC. 

We model the nuclei as uniform spheres and assume a 
transverse profile for the initial temperature as T(r) — 
ro[2(l-rVi?2^)]i/4 where Ri_ = 1.2A^/^ fm. We use the 
same profile R{r) = 2(1 — / R\) for the jet production 
in @ while performing the space-time integration d'^x = 
rdr rdr drjdcj). The limits of the r-integration are [to,t/] 
where Tf is fixed through the relation T^t = const, when 
the temperature reaches Tf = 160 MeV. 

We also give results for the competing photon pro- 
cesses, the direct production in primary hard annihi- 
lation or Compton processes and the production via 
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FIG. 1: Spectrum dN/d^p±dy of photons at y = for cen- 
tral collision of gold nuclei at V^ivjv = 200 GeV at RHIC. We 
show the photons from jets interacting with the medium (solid 
line) , direct hard photons (long dashed) , bremsstrahlung pho- 
tons (short dashed) and thermal photons (dotted). 
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FIG. 2: The same as Fig.Qfor central collision of lead nuclei 
at y^%^ = 5.5 TeV at LHC. 



bremsstrahlung from a produced jet parton. For both 
cases the leading order results can again be found in j^. 
No i^T-factor is used for direct photon production. Our 
results here obey the scaling for single photons suggested 
m from empirical considerations. 

In Fig. n] we plot the results for thermal photons, direct 
photons due to primary processes, bremsstrahlung pho- 
tons and the photons coming from jets passing though the 
QGP in central collision of gold nuclei at RHIC energies. 
The corresponding results for LHC energies are shown in 
Fig. El We see that the quark jets passing through the 
QGP give rise to a large yield of high energy photons. 
This contribution should be absent in pp collisions. For 
RHIC this contribution is the dominant source of photons 
up to p_L « 6 GeV. Note that the jet-to-photon conver- 
sion falls more rapidly with p± than the direct photon 
yiield, similar to a higher twist correction. 

A suppression of the bremsstrahlung contribution due 
to the multiple scattering suffered by the fragmenting 
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partons |27| will further enhance the importance of the 
jet-photon conversion process. The effects of a possible 
chemical non-equilibrium may be estimated as follows. A 
reduced quark abundance will enter at least linearly into 
the thermal rate but can reduce the jet-photon conversion 
by at most a factor of two, as the Compton contribution 
will remain unaffected. Therefore the thermal rate will be 
more strongly suppressed than the jet-photon conversion. 
Similarly, the jet-photon conversion is only moderately 
sensitive to the initial conditions and the transit time of 
the jet in the medium. We have found that increasing 
To to 1 fm/c, keeping TqTq fixed, reduces our results for 
RHIC only by about 40 percent. 

One may ask how our results will be modified by the 
energy loss experienced by the quark jets passing through 
the QGP. This effect will be relatively small for two rea- 
sons. First, we are concerned here with quark jets which 
lose energy at less than half the rate as gluons, which 
form the dominant fraction of jet events. Second, as the 
energy loss scales as the square of the distance traveled 
through the medium 0, the average energy loss at the 
moment of jet-photon conversion is only one-third of the 
total normal energy loss. This argument implies that, 
while the high-p^ hadron spectra encode information 
about the energy loss in the plasma, the photons from 
medium induced emission, on average, carry information 
about jets at an earlier stage. Obviously, the compari- 
son of the high-pj^ hadron and direct photon spectra can 
shed light on the evolution of the jet spectrum. We be- 
lieve that this could be a useful method for the analysis 
of jet quenching and medium effects on hard particles in 
the QGP. More details of the calculation together with 
a careful analysis of the assumptions made here will be 
presented in a forthcoming paper. 

In summary, we have discussed a novel mechanism 
for the production of high energy photons in relativistic 
heavy ion collisions. The photons arise from the passage 
of quark jets through the medium, when a quark is con- 
verted into a photon due to Compton scattering from a 
gluon or annihilation by an antiquark. Our calculations 
indicate that this mechanism is the leading source of di- 
rectly produced photons at RHIC in the region p± < 6 
GeV/c. The spectrum of the resulting photons reflects 
the initial distribution of the hard scattered quarks and 
is sensitive to the initial conditions and the traversal time 
of the jet in the plasma. 
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